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ABSTRACT 

We present the first ultraviolet (UV) and multi-epoch optical spectroscopy of 30 Dor 016, a massive 
02-type star on the periphery of 30 Doradus in the Large Magellanic Cloud. The UV data were 
obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope as part of the Servicing 
Mission Observatory Verification program after Servicing Mission 4, and reveal #016 to have one of 
the fastest stellar winds known. From analysis of the C IV AA1548-51 doublet we find a terminal 
velocity, Voo = 3450 ± 50kms _1 . Optical spectroscopy is from the VLT-FLAMES Tarantula Survey, 
from which we rule out a massive companion (with 2 d < P < 1 yr) to a confidence of 98%. The radial 
velocity of #016 is offset from the systemic value by — 85kms _1 , suggesting that the star has traveled 
the 120 pc from the core of 30 Doradus as a runaway, ejected via dynamical interactions. 
Subject headings: open clusters and associations: individual (30 Doradus) — stars: early-type — stars: 
fundamental parameters — stars: mass- loss 



1. INTRODUCTION 

30 Doradus in the Large Magellanic Cloud (LMC) is 
the richest Hll region in the Local Group, providing an 
excellent template with which to study regions of intense 
star formation, and both stellar and cluster evolution. 
It harbors a significant fraction of the most massive and 
luminous stars known, with a rich population of the ear- 
liest O-type stars (e.g., Melnick 1985; Walborn & Blades 
1997), particularly in R136, the dense cluster at its core 
(Massey & Hunter 1998). 

The 03 spectral class was introduced by Walborn 
(1971) to accommodate stars in which He I A4471 was ab- 
sent in moderate-resolution photographic spectra, com- 
pared to the very weak absorption seen in 04-type spec- 
tra. The classification scheme was extended further by 
Walborn et al. (2002) to include the new types of 02 
and 03.5 to delineate the behavior of the N III, N IV, 
and N V features seen in the earliest types. Even in the 
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age of large multi-object surveys only a few tens of stars 
are known with 02-03.5 types. Although very rare, 
their influence is far-reaching as they are expected to 
evolve rapidly into nitrogen rich Wolf-Rayet stars (WN 
types), plausible progenitors of supernovae and, poten- 
tially, gamma-ray bursts (Smartt 2009). 

Observations with the 2-degree Field (2dF) instrument 
at the Anglo- Australian Telescope revealed a new 02- 
type star on the western fringes of 30 Doradus (Figure [T]), 
with a radial velocity of ~85kms _1 lower than the sys- 
temic velocity of nearby massive stars (e.g., Bosch, Ter- 
levich & Terlevich 2009). New multi-epoch spectroscopy 
of this star, 30 Dor 016 in the VLT-FLAMES Taran- 
tula Survey (Evans et al. 2010), now enables us to rule 
out the presence of a close massive companion to a high 
level of confidence, suggesting that the star might have 
been ejected from the denser central region. Here, we 
combine these observations with new ultraviolet (UV) 
spectroscopy of #016, some of the first data taken with 
the Cosmic Origins Spectrograph (COS) on the Hubble 
Space Telescope (HST), which reveal the star to have one 
of the highest wind terminal velocities seen to date in any 
massive star. 

2. OBSERVATIONS 

2.1. UV spectroscopy 

30 Dor 016 (a = 05 h 37 m 08*88, 5 = -69° 07' 20 / /36, 
J2000) was observed with HST-COS as part of the Ser- 
vicing Mission Observatory Verification (SMOV) pro- 
gram in 2009 July using the G130M and G160M gratings 
in the far-UV channel. The data were obtained as part of 
a focus-check run early in the SMOV phase of COS oper- 
ations so we do not have full wavelength coverage as the 
gratings were not moved. Out-of- focus exposures (<25% 
of the total) were excluded when co-adding the spectra, 
with final exposure times of 5671s (G130M) and 4511s 
(G160M). Estimates of the signal-to-noise ratio (S/N) 
in the combined spectrum are complicated by the large 
number of metallic absorption lines in the spectra of hot 
stars, particularly while the instrument is still undergo- 
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Table 1 

Observational epochs and differential radial velocities 



#016 





7 + 



Figure 1. Digital Sky Survey "IR" (~7-band) image showing the 
location of 30 Dor 016 (encircled) relative to the 30 Dor complex. 



ing full characterization. Here, we use selected contin- 
uum regions that are relatively free of lines to estimate 
the noise, yielding S/N >40. 

Preliminary characterization of the COS on-orbit per- 
formance has found broad non- Gaussian wings in the 
line-spread function of the spectrograph owing to wave- 
front errors from the primary and secondary mirrors of 
the telescope (Ghavamian et al. 2009). This reduces 
the detection threshold for weak absorption features and 
modifies the profiles of saturated lines. To account for 
these effects, our models of the C IV line in Section [3] 
were convolved by the line-spread functions determined 
by the COS team. 

The COS spectrum of #016 is shown in Figure [2j to- 
gether with an HST Space Telescope Imaging Spectro- 
graph (STIS) E140M spectrum from program 9434 (P.I. 
Lauroesch) of HDE 269810 (Sk-67° 211), classified as 
02 III(f*) by Walborn et al. (2002). They are close 
twins, but with subtle differences — the most striking 
being the morphology of the C IV absorption and the 
P Cygni emission at He II A1640. Indeed, the C IV dou- 
blet in #016 is saturated to a large velocity, indicative 
of a very fast stellar wind. 

2.2. Optical data 

Optical spectroscopy of #016 was obtained as part of 
two programs with 2dF, before and after the AAOmega 
upgrade (Sharp et al., 2006). Three 2dF- AAOmega spec- 
tra were obtained on 2006 February 22-23 (P.I. JVL), 
with a 2dF spectrum obtained on 2004 December 8 (P.I. 
IDH), as summarized in Table [I] 

The star was also observed as part of the Tarantula 
Survey (full details to be published elsewhere). Obser- 
vations of #016 with three of the standard FLAMES 
settings (LR02, LR03, and HR15N) were obtained over 
the period 2008 December 17-22. Pairs of exposures 
were taken in each observing block (OB), with three 
OBs observed at the LR02 and LR03 settings, and two 
at HR15N. Individual exposure times were 1815s (LR02, 
LR03) and 2265s (HR15N). The spectral coverage and 



Instrument 


FWHM 


Epoch 


MJD 




and Setting 


[A] 






[kms- 1 ] 




0.61 


1 


r /|o-| 7 OOQ 

54817. 116 


0.0 


[AAoy0U-4004J 




Z 


04ol I .ZD I 


1 7 






3 


04:0ZiZi.U00 


— 1.8 






4 


54860.105 


-1.7 






5 


54890.041 


-1.7 






6 


55112.361 


-2.0 


FLAMES-LR03 


0.56 


1 


54818.243 


0.5 


[AA4499-5071] 




2 


54818.287 


-0.4 






3 


54818.330 


-2.8 


FLAMES-HR15N 


0.41 


1 


54818.076 




[AA6442-6817] 




2 


54818.129 




2dF-AAOmega 


1.0 




53788.442 


8.8 


[1700B] 


1.0 




53788.490 


-3.0 


[1500V] 


1.25 




53789.430 


-3.6 


2dF [1200BJ 


2.75 




53347.564 





Note. — Av r results from cross-correlation with the first epoch 
LR02 observations. 

resolution, as defined by the full width at half- maximum 
(FWHM) of the arc calibration lines, is given in Table [TJ 
The survey features repeat LR02 observations to detect 
radial velocity variables; for #016 these were obtained on 
2009 January 29, February 28, and October 8 (epochs 4, 
5, and 6 in Table [l] respectively). The modified Julian 
Dates (MJD) for each epoch are listed in Table [I] Each 
co-added pair has an S/N > 100. 



2.2.1. Spectral classification 

Informed by the apparent absen ce of detectable radial- 
velocity variations (Section 2.2.2[ ), the individual LR02 
and LR03 spectra were co- added, then merged in the 
overlap region, yielding an S/N > 350. 

The AA3950-4750A region of #016 is shown in Fig- 
ure [2j compared with spectra of HDE 269810 and 
HD 93129A from Walborn et al. (2002). The most no- 
table feature in the spectrum of #016 is the strong N IV 
A4058 emission combined with an absence of N III A4640 
emission, suggesting that the star is of the very earliest 
type in the existing classification framework. Also note 
the O IV A4632, Si IV A4654, and C IV A4658 emission, 
and the He II A4686 P Cygni profile. Accurate rectifica- 
tion of the A4686 and Ka regions is a notorious problem 
in observations of O-type stars, but we see broadened 
wings to the He II line in both manual and scripted rec- 
tifications. 

As in the UV, the spectrum of #016 is similar to that 
of HDE 269810. While the He II A4686 profile in #016 
suggests a somewhat more luminous star compared to 
HDE 269810, it is not as remarkable as the emission seen 
in HD93129A (02 IP; Walborn et al. 2002), leading to 
our intermediate classification of 02 III-IP. 



2.2.2. Stellar radial velocities 

There are relatively few optical lines for radial veloc- 
ity analysis in such an early-type spectrum. The Balmer 
hydrogen lines and He II A4686 are strongly affected by 
the wind, so do not provide a true diagnostic of the stel- 
lar radial velocity. Overplotting of the data from each 
epoch (including those from 2dF) revealed no significant 
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Figure 2. Upper panel: HST-COS spectrum of 30 Dor 016 and an HST-STIS spectrum of its near twin, HDE 269810 (both binned to 
0.25 A). The wind features identified are N V AA1239-43, O IV AA1339-44, O V A1371, C IV AA1548-51, He II A1640, and N IV A1718. Lower 
panel: FLAMES spectrum of #016, illustrating its place in the luminosity sequence at type 02 compared to HD93129A and HDE 269810. 
Emission lines identified in #016 are N IV A4058; Si IV AA4089, 4116, 4654; O IV A4632; C IV A4658. Absorption lines identified in 
HDE 269810 are He II AA4026, 4200, 4542, 4686; N V AA4604, 4620. The cmfgen model for #016 is overplotted in red. 



velocity shifts nor broadening/ asymmetries in the lines, 
suggesting that a high-mass companion is not present 
and that #016 is a prime candidate as an ejected, mas- 
sive runaway from R136. 

The XCORR routine in the DIPSO spectral analysis pack- 
age (Howarth et al. 2004) was used to cross-correlate 
the co-added pairs of LR02 spectra with the first epoch 
of LR02 observations. The velocity offsets, Av r (see Ta- 
ble [IJ, are found from fits to the cross-correlation func- 
tions (CCFs) with Gaussian profiles generated using the 
DIPSO emission-line fitting (ELF) routine. Although the 



CCFs are not Gaussian, the ELF fits were restricted to 
the range —100 < Av r < 100 kms -1 , yielding robust esti- 
mates with fitting errors of 3-4 kms -1 (i.e., one tenth of a 
resolution element). From ELF fits to the He II AA4026, 
4200, 4542 absorption, and Si IV AA4098, 4116 emission 
lines, the mean heliocentric radial velocity for the first 
LR02 observations was v T = 191.8 ±1.4 kms -1 . 

Using the same methods we also cross-correlated the 
relevant calibration arcs to investigate the wavelength 
stability of FLAMES. The instrument is remarkably sta- 
ble over all of our observations, with a largest offset of 
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Av r = — 0.3 kms -1 , found between the arcs for epochs 1 
and 6 (i.e., separated by 10 months). 

Table [T] also lists cross-correlation results for the re- 
gions of the LR03 and 2dF-AAOmega spectra which 
overlap with the LR02 data; good agreement is found 
within the velocity resolution of the respective data. 

2.2.3. Limits on a binary companion 

To investigate the confidence to which we can rule out 
a close massive companion to #016, we used the method 
of Sana, Gosset & Evans (2009). In these calculations 
we adopted a 75 M primary, a bi-uniform period (P) 
distribution (with 50% of the simulated systems with 
0.3 < log P(d) < 1 and with 50% with 1 < log P(d) < 3.5) , 
uniform distributions of eccentricities and mass ratios 
(between and 0.9, and 0.1 and 1.0, respectively), a 
random orientation in space, and a random periastron 
passage time. 

Our observations allow us to rule out an O- or 
early B-type companion (with peak-to-peak variations 
> 20 kms -1 ) with a period in the range 2d < P < lyr at 
the 98.1% level. These results are only weakly depen- 
dent on the primary mass, with a variation of ±1% for a 
primary in the range of 50-90 M . 

3. TERMINAL WIND VELOCITY 

The Sobolev with exact integration (SEI) method was 
used to calculate UV line profiles to estimate v^, the 
terminal wind velocity of #016. We used the CCP7 
SEI cod a^ which is based on the method from Lamers, 
Cerruti-Sola & Perinotto (1987). The observed C IV 
AA1 548-51 doublet was blueshifted by 190 kms -1 (Sec- 
tion 2.2.2), then compared to the calculated profiles. 



The code adopts a standard /3-law to describe the 
acceleration of the stellar wind, and includes a turbu- 
lence parameter, v t (typically expressed as a fraction of 
Vqq), to account for stochastic deviations from the ve- 
locity prescription. Line profiles were calculated over 
a range of values of Voo, /?, and u t , then con volved by 
the line-spread function noted in Section |2.1| The best 
X-squared fit (excluding the interstellar absorption fea- 
tures and the emission peak, where the SEI method 
is known to underestimate the flux by about 10%) 
is obtained with: ^ = 3450 ±50 kms -1 , f3 = 0.75 ±0.1, 
v t = 0.02 ± 0.01 Voq. The optical depth law parameters 
in the code (cti and o^) were set to 1. For saturated 
lines these can be varied to obtain subtly different fits to 
the redward part of the absorption component but with 
no significant impact on the derived velocities. When v t 
is relatively small, the value of Voo obtained from fitting 
the line profile effectively corresponds to the blueward 
extent (i.e. ^biack) of the saturated P Cygni profile. 

This is one of the highest wind terminal veloci- 
ties measured directly for an O-type star, behind only 
HDE269810 (v^ = 3750 kms" 1 ; Walborn et al. 1995) 
and R136a-60£p] (03 IP, = 3640 kms" 1 ; Prinja & 
Crowther 1998^ 



11 http://ccp7.dur.ac.uk/library.html 

12 Star 36 from Massey & Hunter (1998). 



4. QUANTITATIVE ANALYSIS 

Physical properties for a subset of 02-type stars were 
presented by Walborn et al. (2004). We employed simi- 
lar methods to inform our discussion of #016, using CM- 
fgen (Hillier & Miller 1998) which solves the radiative 
transfer equation in the co-moving frame, under the addi- 
tional constraint of statistical equilibrium. The temper- 
ature structure follows from the assumption of radiative 
equilibrium. CMFGEN does not currently solve the mo- 
mentum equation, so a density or velocity structure is re- 
quired. The velocity for the supersonic part of the wind is 
parameterized with a classical /3-type law (f3 = 0.75, from 
Section [3| , which is connected to a hydrostatic density 
structure at depth, such that the velocity and its gradi- 
ent match at the interface. The subsonic density struc- 
ture is set using a line-blanketed, plane-parallel tlusty 
model (v.200; Lanz & Hubeny 2003) with log # = 3.75 
(preferred over log # = 4.0). The atomic model is similar 
to that adopted by Walborn et al. (2004), including ions 
from H, He, C, N, O, Ne, Si, P, S, Ar, Fe and Ni. 

We assumed a depth-independent Doppler profile for 
all lines when solving for the atmospheric structure in 
the co-moving frame, and incoherent electron scatter- 
ing and Stark broadening for hydrogen and helium lines 
was included. A uniform turbulence of 50 kms -1 was 
adopted in the calculation of the emergent spectrum in 
the observer's frame. Using the same methods as Wal- 
born et al. (2010), a line-broadening parameter (y sini) 
of 88±17kms _1 was found from the optical spectra. 
Clumping was incorporated using a volume filling factor, 
/, as described by Hillier et al. (2003), with a typical 
value of / = 0.1 resulting in a reduction in mass-loss rate 
by a factor of ~ 3. 

4.1. Spectroscopic results 

Based upon the observed UV to near-IR spectral en- 
ergy distribution and the Ha and He II AA1640, 4686 
wind features (with = 3450 kms -1 ), we estimated the 
stellar temperature, luminosity and mass-loss rate of 
#016 simultaneously. 

For an adopted LMC distance of 49 kpc (distance mod- 
ulus = 18.45) the UV spectrophotometry and opticaj^j 
and near-IR photometr^] of #016 can be reproduced 
with interstellar extinctions of E(B — V)=0.07 (fore- 
ground) and 0.27 (LMC), using standard Galactic and 
LMC (field) extinction laws (Seaton 1979; Howarth 1983) 
and a stellar luminosity of L* = 1.2 x 10 6 L , as shown in 
Figure [3J The relatively large local extinction is sup- 
ported by the strong A4428 diffuse interstellar band, 
which is at a comparable radial velocity to #016. 

In the absence of measurable He I A4471 absorption, a 
combination of optical nitrogen lines (N IV A4058, N V 
AA4603-20) and UV oxygen lines (O IV AA1339-44, O V 
A 1371) was used to derive an estimated effective tem- 
perature (T e ff) of 50 kK for #016. This is a compromise 
between the available diagnostics — the optical N V lines 
favor a temperature some 5%-10% higher, while the UV 
O IV suggests 10% lower — but the strong N IV A1718 
P Cygni profile and negligible S V A 1501 line both firmly 

13 U= 12.62, B = 13.53, V = 13.49, # = 13.43 (Massey 2002). 

14 J = 13.38, H = 13.35, K s = 13.36 (Two Micron All Sky Survey, 
Skrutskie et al. 2006) 
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support T e ff = 50 kK. Given the adopted wind parame- 
ters, a slightly greater broadening profile was preferred 
(v sini ~150kms _1 ) in the final model than determined 
above. 

A range of nitrogen abundances were considered, 
adopting Xn = 0.08% by mass, which is substantially 
greater than the interstellar abundance in the LMC 
(0.016%; Russell & Dopita 1990). For carbon and oxygen 
we adopt Xq =0.08% and Xq =0.2% by mass, slightly 
below the results from Dufour (1984) and Russell & 
Dopita (1990) for the interstellar medium (0.095% and 
0.35%, respectively), although the UV oxygen diagnos- 
tics are relatively insensitive to abundance changes. 

The adopted model is shown in Figures [2] and [3j in- 
cluding allowance for an LMC neutral hydrogen column 
density of 10 21,6 cm -2 . In general, the fits are satisfac- 
tory, with a few exceptions. The UV P Cygni absorption 
troughs of O V A1371 and N IV A1718 extend too far blue- 
ward, although the majority of the "iron forest" (Fe V- 
Vl) features are well matched. As with the temperature, 
the adopted mass- loss rate (M = 3 x 1O _6 M yr _1 ) is 
a compromise between the available diagnostics (He II 
AA1640, 4686 and Ha). Weak metallic lines in the opti- 
cal are, in general, reproduced satisfactorily, e.g., Si IV 
AA4089-4116, O IV A4632, and the blend of C IV A4658 
and O IV A4654, 4663. A more detailed treatment of 
#016 (including the role of clumping, line-broadening, 
and the wind acceleration and turbulence) is beyond the 
scope of this Letter and will be presented elsewhere. 

A summary of physical and wind parameters of 
#016 is provided in Table [2j compared with those for 
HDE 269810 from Walborn et al. (2004). We have also 
compared our results with LMC-metallicity, non-rotating 
evolutionary predictions (R. Hirschi, 2010, private com- 
munication) calculated with the mass-loss recipe of Vink, 
de Koter & Lamers (2001). This approach is supported 
by our derived mass-loss rate, which agrees with the- 
oretical predictions to within ~0.05dex for #016. We 
obtain the best agreement for a stellar mass of ~90 M 
at an age of ~1 Myr, corresponding to log g ~ 4.0. In con- 
trast, the spectroscopic gravity suggests a considerably 
lower mass of ~ 50 M . Using the evolutionary mass to 
calculate the escape velocity, v esc , our measured termi- 
nal velocity is in excellent agreement with the expected 
value using the empirical relation ~2.6v esc (Lamers, 
Snow & Lindholm 1995), as well as new predictions for 
the velocity structure of massive star winds (L. Muijres 
et al. 2010, in preparation). 



5. 30 DOR 016 AS A MASSIVE RUNAWAY 

30 Dor 016 is located at the end of an extended rar- 
efied filament of nebular gas, at a radial distance of ~8f25 
(120 pc in projection) west-southwest of R136, and 4f75 
(70 pc) from the less dense cluster NGC 2060. The [O III] 
emission superimposed on the stellar spectrum has a 
peak velocity of ~275kms _1 , with some asymmetry to- 
ward slightly lower velocities but with no evidence for a 
second component matching the velocity of #016. This 
suggests #016 did not form locally and is a bona fide 
runaway star. 

Further analysis of the stellar and gas dynamics in this 
region is warranted to complete our picture of the origins 
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Note. — Mass estimates are obtained from non-rotating evo- 
lutionary models at LMC metallicity (R. Hirschi, private commu- 
nication). Metal abundances are reliable to within 50% (N), or a 
factor of 2 (C and O). 



of #016, e.g., van Loon & Zijlstra (2001) found molec- 
ular gas to the southeast of R136 at a comparable ra- 
dial velocity to the star. At present we lack proper mo- 
tions at useful precision to constrain its transverse veloc- 
ity: n a = 1.2 ± 5.9, fji 6 = 1.5 ± 5.9 mas yr" 1 from UCAC2; 
/i a = 2.8±3.8, /i 5 = -4.3±3.8 mas yr" 1 from UCAC3 
(and flagged as uncertain; Zacharias et al. 2010). HST 
imaging of this region would provide much improved con- 
straints (e.g., the precision achieved by Kallivayalil et al. 
2006), and in the longer term, the Gaia mission should 
deliver transverse velocities in the LMC to a few kms -1 . 

There are two other 02-type stars (Sk— 68° 137 and 
BI 253) to the north and northwest of 30 Dor, which were 
proposed as runaways by Walborn et al. (2002). BI253 
has also been observed within the Tarantula Survey, but 
from preliminary inspection of its spectra, its radial ve- 
locity appears consistent with this part of the LMC. The 
case of #016 is also reminiscent of the discovery of Nll- 
026 (02.5 III(f*)), suggested as a runaway by Evans et 
al. (2006). 

Runaway stars are thought to result from either dy- 
namical interaction in massive dense clusters, or via a 
kick from a supernova explosion in a binary system, with 
the more massive star exploding first (see, e.g., Gvara- 
madze, Gualandris & Portegies Zwart 2009, and refer- 
ences therein). It is generally accepted that R136 is suf- 
ficiently young (1-2 Myr) that its most massive stars have 
yet to explode as supernovae. This implies that, if from 
R136, #016 must have been ejected through dynamical 
interaction, one of the clearest cases to date in support of 
this mechanism. This is vitally important as dynamical 
interactions in massive clusters are thought to be a pos- 
sible mechanism for producing stellar mergers and very 
massive stars (with masses > 140 M ) which might sub- 
sequently end their lives as pair-instability supernovae, of 
broader relevance to the early ages of the universe when 
such massive stars are thought to be common (e.g., Heger 
et al. 2003; Gal- Yam et al. 2009). 
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Figure 3. CMFGEN fits (red) to the flux-calibrated UV and photometric observations (black lines/squares), used to determine the interstellar 
extinction toward #016. Insets are expanded plots of the fits to the HST-COS spectra. 
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